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A NOTE ON ALGEBRAIC RICCATI EQUATIONS ASSOCIATED 
WITH REDUCIBLE SINGULAR M-MATRICES 

DI LU AND CHUN-HUA GUO 


Abstract. We prove a conjecture about the minimal nonnegative solutions 
of algebraic Riccati equations associated with reducible singular M-matrices. 
The result enhances our understanding of the behaviour of doubling algorithms 
for finding the minimal nonnegative solutions. 


1. Introduction 
For the algebraic Riccati equation 


( 1 ) 


XCX - XD-AX + B = 0, 


where A, R, C, D are real matrices of sizes m x m,m x n,n x m,n x n, respectively, 
a systematic study was done in [1] when 


( 2 ) 


K = 


D -C 
-B A 


is a nonsingular M-matrix or an irreducible singular M-matrix. The study was re¬ 
cently extended in [2] to reducible singular M-matrices under suitable assumptions. 

A real square matrix A is called a Z-matrix if all its off-diagonal elements are 
nonpositive, so any Z-matrix A can be written as si — B with R > 0. A Z-matrix A 
is called an M-matrix if s > p{B), where p{-) is the spectral radius; it is a singular 
M-matrix if s = p{B) and a nonsingular M-matrix if s > p{B). 

Some regularity assumption is needed to guarantee the existence of a solution of 
the equation (HD associated with the M-matrix K. An M-matrix A is said to be 
regular if Av > 0 for some u > 0. 

The following result is proved in [5]. 


Theorem 1. Suppose the matrix K in ^ is a regular M-matrix. Then (HD has a 
minimal nonnegative solution $ and D — is a regular M-matrix, and the dual 
equation 


(3) YBY -YA-DY + C = Q, 

has a minimal nonnegative solution and A—Rd' is a regular M-matrix. Moreover, 

Im — and are both regular M-matrices. 
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Associated with the matrix K in © is the matrix 


( 4 ) 


H = 


■ In 

0 

K = 

■ D 

-C ' 

0 

-^m 

B 

-A 


When the matrix AT in ([2|) is regular singular M-matrix, the following assumption 
has been introduced in [5]. 

Assumption 1. The matrix H in o has only one linearly independent eigenvector 
corresponding to the zero eigenvalue of multiplicity r > 1. 

It is known [T] that the assumption is satisfied with r = 2 when K is an irre¬ 
ducible singular M-matrix. 


Under Assumption [TJ there are nonnegative nonzero vectors 
where ui,vi G K" and U 2 ,V 2 G K™, such that 


Ul 

U2 


and 


Vl 

V2 


K 


Vl 

V2 


= 0, [ul ul]K = 0. 


They are each unique up to a scalar multiple [5]. 

The purpose of this note is to provide an afSrmative answer to a conjecture in 
[5], regarding the matrices Im — 4>'I' and /„ — ^'d). 

2. The result 

The following result was conjectured to be true in [2] , and was proved under the 
restrictive assumption that at least one of $ and is positive. 

Theorem 2. Let K he a regular singular M-matrix with Assumption{J[ Ifujvi ^ 
u^V 2 , then Im — 4>'I' and are nonsingular M-matrices. 

Proof. By Theorem[Tl Im — 4>'I' and are both M-matrices. So we just need 

to show they are nonsingular when ujvi u"^V 2 . Since is nonsingular if 

and only if /„ — $'1' is nonsingular, we only need to show Im — is nonsingular. 


In view of 

'In O' 


■ In 



'In ^ ' 


. -4> Im _ 


$ 

Im 


_ 0 

we need to show that the matrix 








■ In 







$ 




' D -C ' 


'In 4 * ' 


'In 4 ' ' 


' R o' 

B -A 


$ Im 


. ^ Im ^ 


_ 0 -s _ 


is nonsingular. Since $ and are solutions of © and respectively, it is easily 
verified that [2] 

( 5 ) 

where R = D — (7$ and S = A — B'it are M-matrices. Therefore, the eigenvalues 
of R and S are all in the closed right half plane, with 0 being the only possible 
eigenvalue on the imaginary axis. When ufvi u'^V 2 , we know from [2] that one 
of the matrices R and S is singular and the other is nonsingular. It follows that 
the matrices R and —S have no eigenvalues in common. 

Let 

' In 4 - 
$ Im 


W = Ker 
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We need to show W = {0}. 

For any x € W, post-multiplying ([5|) by a; shows that Tx G W, where T is the 
linear transformation from to defined by 


yi 


■ R 

0 


yi 

. 2^2 . 


0 

-S' 


. 2/2 _ 


where yi G C" and y2 G C™. 
transformation T. Suppose W 


Thus W is an invariant subspace of the linear 


^ {o: 


Then we have 0 ^ w G W such that 


Tw = Xw for some A G C. Write w = 


Wi 

W2 


where wi G C" and W2 G C"*. We 


then have Rwi = Xwi and {—S)w 2 = Xw 2 - Since R and —S have no eigenvalues in 
common, one of wi and W 2 must be a zero vector. It then follows from 


'In ■ 


Wi 



W2 


that wi and W 2 are both zero vectors. The contradiction shows that W = {0}. □ 


It has been explained in [2] that when K is a regular singular M-matrix with 
Assumption [T] and ujvi ^ u^V 2 : doubling algorithms (such as SDA in [4] and 
ADDA in 0) can be used to find the minimal nonnegative solutions $ and 'I' 
simultaneously. But before Theorem [2] is proved in this note, there were a few 
subtle issues associated with the doubling algorithms. We now have the following 
modification of [51 Theorem 15] about the ADDA, which uses two parameters a 
and p. The ADDA is reduced to the SDA when a = p. 

Theorem 3. Let K be a regular singular M-matrix with Assumption[J\and ufvi ^ 
U 2 V 2 - Assume that a > maxau > 0 and (3 > maxda > 0. Then the ADDA is 
well defined with I — GkHk and I — H^Gk being nonsingular M-matrices for each 
k > 0. Moreover, Eq < 0,Fo < 0, Ek > 0, Ffc > 0, 0 < Elk-i < Hk < 

0 < Gk-i < Gfc < Ip for all k > 1, and 

limsup ^'{/WHk - $11 < r(a,/3), limsup - 4'|| < r(a,/3), 

k—foo k—foc 

where r{a, fi) = p ((i? -f aI)~^{R — fil)) ■ p (^{S -f — a/)) < 1 with R = 

D-G$,S' = A-ST. 

Since we have now proved that J—and J—Tih are nonsingular M-matrices we 
can use the approach in [3] to prove that the ADDA is well defined with I — GkHk 
and I — HkGk being nonsingular M-matrices for each k > 0 even when a = max an 
and (3 = maxdn. That r{a,j3) < 1 in Theorem|3]is already known in |5]. Thus, Elk 
converges to $ quadratically, and Gk converges to quadratically. 

By 0 Theorem 2.3], the parameters a = max an and (3 = maxdn minimize 
r(a, (3) among all parameters a > maxa^i and ft > maxdn- Therefore, we should 
normally use the optimal values a = maxa^i and (3 = maxdn for the ADDA. Before 
Theorem [5] is proved, we avoid using a = maxo^i and (3 = maxdn, to ensure that 
the ADDA is well defined. 

The matrices (/ — GkHk)~^ and (/ — HkGk)~^ appear in the ADDA. With the 
proof of Theorem|51 we now know that, in Theorem|31 the matrices I — GkHk and 
/ — HkGk will not converge to singular matrices. This is of course favorable for the 
ADDA. 
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